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on short range electron-lattice and spin-orbital correlations we supplement the 
modelling with numerical simulations. 
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1. Introduction 

Mixed-valence manganese oxides of perovskite structure have been studied for more 
than fifty years [1,2]. However, many of their peculiar properties are still rather 
poorly understood. About a decade ago the observation of the so-called colossal 
magneto-resistance effect (CMR) [3-5] moved these materials again into the focus of 
intense research activity [6-8]. It turned out soon that the complex electronic and 
magnetic properties of the manganites as well as their rich phase diagram depend 
on a close interplay of almost all degrees of freedom known in solid state physics, 
namely itinerant charges, localised spins, electronic orbitals, and lattice vibrations. 
To separate important and less important degrees of freedom and to work out the 
essential low energy physics are therefore two particularly crucial factors for any 
microscopic description of these interesting compounds. In the present article we 
give a detailed introduction to the electronic structure of manganites and derive a 
microscopic model, which includes the dynamics of charge, spin, orbital and lattice 
degrees of freedom on a quantum mechanical level. We also discuss potential further 
simplifications and complement the work with a review of previous numerical studies 
of the derived model. 



Figure 1. Left: Idealised perovskite structure (space group Pm3m). Right: Real 
structure of LaMnC>3 at T = 1.5 K (orthorhombic, space group Pnma) [9]. 



Table 1. The four symmetry operations that generate the cubic point group Oh, 
the corresponding notation as well as the action on the position vector. 
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Symbol 


Coordinate transformation 


inversion 
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{x,y,z) 


-» {-x, -y, -z) 


4-fold rotation about x-axis 


Cf 


{x,y,z) 


■* (x,-z,y) 


4-fold rotation about z-axis 


Cf 


{x,y,z) 


■* {-y,x,z) 


3-fold rotation about diagonal 


cl 


{x,y,z) 


•+ (z,x,y) 



2. Crystal structure and symmetries 

The structural phase diagrams observed for CMR manganites are almost as rich as 
their electronic and magnetic counterparts (for a review see, e.g., Ref. [6]). Consider, 
for example, the undoped parent compound LaMnOa. At low temperature its 
structure is characterised by the orthorhombic space group Pnma [9] , but doping with 
strontium, Lai-^Sr^MnOs, will transform the system to the triclinic space group R3c 
at i k 0.175 [9,10]. On the other hand, similar structural transitions occur with 
increasing temperature [11]. However, nearly all of these crystal arrangements can be 
understood in terms of distortions of the ideal perovskite structure shown in Figure^ 
Here each manganese ion is surrounded by an octahedron of oxygen atoms, whereas 
the large rare earth (R) or alkaline earth atoms (A), used as dopants, occupy the 
centre of the cube formed by the manganese sites. 

For a theoretical description it is thus quite natural to start from the above 
ideal cubic structure, and to account for deviations by including (dynamical) 
lattice distortions, i.e., phonons, into the microscopic model. Since the cubic site 
symmetry of the manganese is particularly important, and since its properties and the 
corresponding notation are frequently used throughout the article, let us first recall 
some of the basic features of the cubic point symmetry group Oh- The group Oh 
consists of 48 symmetry operations which can all be generated by the inversion of 
space and the 3 basic rotations listed in Table Q Since the rotations commute with 
the inversion the group Oh can be represented as a direct product, Oh — O (8> Cj, of 
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Table 2. The five irreducible representations of the cubic group O, the notation 
used for standard basis elements, examples of corresponding functions, and 
transformation properties with respect to the group generators. Note that we 
rotate the functions, not the coordinate system. 
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the inversion group Ci and the group O, which is formed by the rotations only. With 
respect to inversion every function can be decomposed into an even (g) and an odd 
(u) part, which reflects the two irreducible representations of Ci with characters 1 and 
— 1, respectively. Similarly, every irreducible representation of O leads to an even and 
an odd irreducible representation of Oh- Given the five irreducible representations of 
O listed in Table |3 we thus find ten irreducible representations of the full cubic point 
group Oh, denoted by Ai g , Ai u , A-i g , etc. Here we follow the notation commonly used 
in the literature (see, e.g., the textbook of Griffith [12]). 

As we know, for example from the coupling of angular momenta, the product of 
elements belonging to two irreducible representations in general leads to a function 
that belongs to a reducible representation, which can again be decomposed into 
irreducible representations. For spins the corresponding coupling coefficients are 
known as the Clebsch Gordan coefficients, and, of course, for the cubic group analogous 
coefficients were calculated and are listed, e.g., in Ref. [12]. For illustration consider 
the product E ®T\, which leads to a six-dimensional reducible representation that 
can be decomposed into T\ and T2. In terms of the above basis functions we find 

, 1 . y/3 v3 . 1 

x = ——6®x-\ — — e <8> x , £ = — — ® x — — s ® x , 

V = --8®y- — e®y, ^=—6®y--e®y, (1) 

z' = 6® z, (' = s(g> z, 

and obviously the new (dashed) functions fulfil the transformation properties of basis 
functions of T\ and Ti- 

3. Coulomb interaction and local electronic structure 

Considering the end members of the manganite series Ri-^A^MnOs at doping x = 
and 1, respectively, within a simplified ionic description the different constituents are 
assigned the formal valences R 3+ Mn 3+ 03~ and A 2+ Mn 4+ Og~. This corresponds 



.4/ 3 



l2g 




A cf S=1/2 

i 



JT 

x>0 



S=3/2 



Mn 



3+/4+ 



Figure 2. Left: The local electronic structure of Mn 3d electrons in a cubic 
crystal field. Right: Angular structure of the e g and t2 g single electron wave 
functions. 



to always completely filled oxygen p-bands and, with increasing x, to hole doping 
of the manganese d-shell, <i 4 — > d 3 . To some degree already the early doping 
dependent measurements of the local magnetic moment in Lai-^Ca^MnOs by 
Wollan and Koehler [2] confirm this general picture, but also recent band structure 
calculations [13, 14] and spectroscopic experiments [15-17] are indicative of manganese 
3d conduction bands. Although these bands are partially subject to hybridisation with 
neighbouring oxygen 2p states, it is thus reasonable and common to consider Mn d 
electrons as the starting point of a microscopic modelling, and proceed along the line 
which is called the weak-field coupling scheme [12]. 

The spherical harmonics Yi. m with I = 2, which describe the angular part of 
the d shell, are all even with respect to inversions. Comparing with Table |3 we thus 
find that these functions can be combined to yield basis functions of the irreducible 
representations E g and T2 g of Oh, 
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Consequently, within the cubic crystal field the five degenerate d levels split into an 
Eg orbital doublet and a T^ g orbital triplet. Taking into account the spatial structure 
of the above wavefunctions and the overlap with neighbouring oxygen electrons, we 
find that the ti g orbitals are lower in energy compared to e g . The magnitude of this 
crystal field splitting, A c f, varies with the occupancy of the levels, typical values for 
the manganites are of the order of 1.5 to 2.5 eV [15]. 

So far we considered only the single-electron picture, but for the manganites 
the on-site Coulomb interaction is strong and thus of particular importance. Before 
describing the many-electron states of manganese ions in detail, let us review some of 
the basic facts for the Coulomb interaction of d electrons. Expanding the Coulomb 
interaction between two charges in terms of spherical harmonics, 

p 2 2 00 / \fc 4 k 

— p — — k 
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with r> = max(ri, r 2 ) and r< = min(ri, r 2 ), we find that the matrix element between 
different single electron functions ip a = R na i a Yi aTna s r7a is given by 

oo 
fe=0 

F a/3lS ={Rn a l a Rn l3 l l3 \—^f\Rn^Rn s l 5 ) 1 (4) 



c fm;i'm' _ \j 2k + 1 I ^ hm - m ' ^J'm' ) • 

Here we already made use of the fact that the Coulomb interaction preserves the z 
component of the angular momentum. However, the special structure of the matrix 
elements cf m . [lm , imposes further restrictions on the expansion index k and suggests 
a close relation to Clebsch Gordan coefficients. Indeed, long ago Racah [18] derived 



c\ m , Vml = {-iy-\j\{2v + i)c lkV ( ^ 
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(5) 



for l + l' + k odd, 

Cikv = < 2(1 + 1' -k)\(l + k-l')l(l' + k-l)\g\ 2 , „ , , 

77"^ 5 TT77 ni9 ; TT^PT, nn for 2g := I + I' + k even , 

(I + I' + k + l)\(g - l)\ 2 (g - 2 (g - k)\ 2 

which leads to the conditions of even l + l' + k and \l — l'\ < k < 1 + 1'. For the present 
case of 3c? electrons the values of k are thus restricted to k = 0, 2, 4, and the Coulomb 
matrix elements reduce to a sum of only three terms, 

2 

K = 

In realistic situations the integrals over the radial parts of the wave functions, 

2 k 

F k = <i?32(ri)i?32(r2)|^|i?32(ri)i?32(r2)) , (7) 

r > 

are usually hard to calculate exactly and are therefore taken as free parameters which 
depend on details of the considered ions and compounds. Instead of the above F k the 
so-called Racah parameters, 

A = F°-F 4 /9, B=— — , C=— , (8) 

' 5 44i 63 ' w 

are more common, because they simplify the notation of the Coulomb matrix elements. 
Since the F k are positive and monotonously decreasing in k [12, 19], the Racah 
parameters are positive as well. In addition, for many transition metals the ratio 
C/B is almost constant and of the order of 4 to 5. Later we make use of this by 
expressing the Coulomb interaction in terms of just two parameters U and J/j. In 
Table we show estimates of A, B and C for manganese obtained from spectroscopic 
data. 

Based on the above introduction we are now in the position to construct many- 
electron cigenstates of the Coulomb interaction for a single ion in a cubic crystal 
field, which define the starting point of our modelling. Neglecting spin-orbit coupling, 
which is usually very small for compounds of the first transition metal series, the 
Coulomb matrix can be decomposed into blocks of given spin and cubic symmetry. 
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Table 3. Different estimates for the Racah parameters of manganese ions (in eV) 
and corresponding references. 



ion 


A 


B 


C 


reference 


Mn 2 + 


6.05 


0.107 


0.477 


[20,21] 




5.43 


0.1190 


0.4122 


[22] 


Mn 3 + 


6.40 


0.120 


0.552 


[20,21] 




5.26 


0.120 


0.552 


[21,23] 


Mn 4 + 


6.58 


0.132 


0.610 


[20,21] 



Especially for the low-energy states it turns out that this symmetrisation diagonalises 
the Coulomb matrix, i.e. we immediately find the ground-states of the n electron 
systems. Taking into account the crystal field splitting, for the Mn 4+ ion id 3 ) the 
ground state is obtained by triply occupying the ti g levels and forming a state of 
maximum spin (Hunds rule). This leads to the spin quartet and orbital singlet A2, 
which in terms of the fermionic creation operators S V(J reads 



and its Coulomb energy is e( 4 A2) = 3A — 15B. For Mn 3+ (d 4 ) the ground state is a 
spin quintet and orbital doublet 5 E with components 



l*> 2 >™> - + \/4|T^^ +)(2+m) ^ c a^ c a |0) ' (10) 



which reflects the freedom of choosing one of the e g levels when adding another electron 
to the d 3 system and requiring strong Hunds rule coupling, i.e. maximal spin. The 
Coulomb energy of the 5 E state is e( 5 E) = 6A — 21B. Note also that in the 9 
component of the many-electron state the e single-electron level is occupied, and vice 
versa, which is caused by the d 3 state belonging to A^. 



4. Perturbation theory in the electronic hopping 

The energy difference between the above ionic ground states and the corresponding 
lowest excitations is given by the crystal field splitting A c f, or by a Coulomb energy 
of at least 8B + 3C. Both energy scales are large compared to the hopping matrix 
clement of e g electrons between two neighbouring manganese sites, which is of the 
order of t w 0.3 to 0.4 eV [24-26]. The hopping of t^g electrons, t w , is even 
smaller. The derivation of an effective low-energy Hamiltonian can thus be based 
on a perturbative expansion in terms of t(~\ . Since the transfer of electrons between 
different manganese sites always proceeds via the completely filled 2p shell of bridging 
oxygen sites, the hopping matrix elements acquire a particular orbital dependence, 
which leads to effective spin-orbital interactions and finally causes all the interesting 
patterns of orbital ordering or disorder observed in the manganites. As Figure 
illustrates for bonds in z-direction, the e orbital by symmetry has no overlap between 
any of the three oxygen p orbitals, whereas for the 9 orbital only the overlap to p z is 
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Figure 3. Illustration of the orbital anisotropy of the electronic hopping via 
bridging oxygen sites. 



finite. Similarly we find that the ti g orbitals have finite overlap at most to one of p x 
or p y , and we can thus summarise all possible hopping processes in z-direction with 
the Hamiltonian 



EM 



i+z,t 



H.c. 



(12) 



At first glance the situation for x and y-bonds appears to be more involved. However, 
using the cubic point symmetry operations, e.g. the diagonal rotation Cf , we can easily 
derive the corresponding hopping matrix elements by simply inserting the rotated 
orbitals into Equation l|12l) . More formally we introduce the operators 

\i d _ (r<d\2 _ (r<d\-i e> _ (jj d y> - 



ci 



Ry - (C; 



3 d ) 2 
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e /e)/2| 
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and obtain for the complete hopping Hamiltonian 
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R& 



tc 



VSfto + ^( C l^cr C 'i+<5^cr + c i, n a C i+S,r]cr) + H.C. 



(14) 



In general i?hopp connects the ionic ground state at one site to a large number 
of Coulomb excited states at the other site. A complete description of the electronic 
subsystem would thus require a huge local Hilbert space, which includes basically 
all ionic excitations. However, these high-energy excitations are irrelevant for most 
of the features we want to model, and we thus restrict the local Hilbert space to 
the ionic ground states derived in Section [21 and account for the excitations in the 
framework of standard degenerate perturbation theory. Then the first order of the 
effective Hamiltonian is given by the matrix elements of -ffhopp between basis states 
of the form 10 and (|10|l on neighbouring sites. These basis states differ in the total 
on-site spin \Si\. To achieve an uniform description it is thus convenient to rewrite 
the spin degrees of freedom in terms of Schwinger bosons [27] , which provide a simple 
means for changing \Si\, 

S t = a i,1 a i,l > S i = K\T a U " 4l a *,lV 2 ' 



' a u a a 



)/2- 



(15) 
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Table 4. Excitations of a single manganese ion that can be reached by adding or 
subtracting an e g or ti g electron to the basis states 5 E and 4 A 2 . Starred Coulomb 
energies in the last column are approximate (see text and Ref. [12]). 

initial state operation final state energy 

5 E[t3( 4 A 2 )e] d 4 / d 5 c\ x 6 Ai[i 3 ( 4 A 2 )e 2 ( 3 A 2 )] 10A - 35B 



4 Ai[f|( 4 A 2 )e 2 ( 3 A 2 )] 10A - 25B + 5C 

4 A 2 [^( 4 A 2 )e 2 ( 1 .4 1 )] 10A-13B + 7C 

4 E[4( 4 A 2 )e 2 ( 1 E)] 10A-21B + 5C** 

cj at 4 Ti[t 4 ( 3 Ti)e] 10A-25B + 6C* 

4 T 2 [t 4 ( 3 Ti)e] 10A-17B + 6C* 

d 4 \d 3 c el 4 A 2 [r|] 3A-15B 

c t2L *Ti[^( s Ti)e] 3A-3B* 

4 T 2 [t 2 ( 3 Ti)e] 3A-15B 

4 A 2 [t 3 ] d 3 /d 4 5 £[t 3 ( 4 A 2 )e] 6A-2LB 

4t 3s [* 2 ( 4 ^2)e] 6A-13B + 4C* 

cj^ 3 Ti[t 4 ] 6A - 15B + 5C* 

d 3 \ d 2 c t | 3 Ti[t 2 ] A - 5B* 



In addition, we can also reduce the number of fcrmionic operators and introduce hole 
operators d\ g and d^l instead. The orbital part of our basis states and corresponding 
projection operators then read 

|0)=4|O>, \s)=4\0), |a 2 )=44|0), 

P e = n g (l-n £ ), Pt = n £ (l-n g ), P a - = n e n g , (16) 

where the vacuum |0) is invariant under cubic symmetry operations, i.e. belongs to 
A\. Calculating the effective hopping matrix element for a z-bond we are mainly 
concerned with the spin part, whereas the orbital part is rather trivial, since only 
9 electrons are allowed to tunnel. Given the spin-2 state \e,2,m) on site i and the 
spin-3/2 state |a2, 3/2, to') on site j the hopping of a 8 electron will transfer a spin-1/2 
from i to j. Projecting the new states onto our basis will lead to I02) at i and |e) at j. 
In addition, the hopping conserves the total bond spin 5j + Sj. Combining all three 
properties of the process, i.e. transfer of a spin-1/2, increase of \SA and lowering of 
I Si |, as well as conservation of S,: + Sj, in terms of Schwinger bosons we can almost 
guess the effective matrix element to be proportional to a'j ^a i ^ + a^a i ^. Indeed, a 
more careful calculation [28] of the double-exchange [29,30] matrix element shows that 
corresponding prefactor is —t/ (25 + 1), where 5 is the amplitude of the shorter of the 
two on-site spins. For a bond in z-direction the first order of the effective Hamiltonian 
is thus given by 

H lf = ~\ ( a i,t4t + a i,l a ll) ''■ "">■ ''.>■■'":■ + Hx - ( 17 ) 
Note that the amplitude of the local spin is restricted by 2|5j| = at ^ + at ^ = 

The second order of perturbation theory in t^) is a bit more involved, since now all 
virtual excitations of the basis states have to be taken into account. In Table 0] we list 



Table 5. The matrix elements (d n d rn \H* opp \d n - 1 d m+1 ) = O a0 S(x)t^ 
together with the corresponding excitation energies Ae. For brevity we denote 
vanishing matrix elements with dots, write x = SiSj and give the squared S(x). 
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all eleven many-particle states that can be constructed by adding or subtracting one 
e g or tig electron to the basis states from Equations 10 and p(J|) . thereby respecting 
cubic and spin symmetry. Since usually there are other configurations belonging to 
the same representations of orbital and spin symmetry, these states are not necessarily 
eigenstates of the on-site Coulomb interaction. However, since these other states have 
no overlap to our basis, taking into account the exact Coulomb eigenstates will only 
slightly modify the energy denominator of the second order terms, but not the general 
structure of the matrix elements. It is thus sufficient to consider the expectation value 
of the Coulomb energy in the above configurations, and the corresponding approximate 
energies are marked with a star in Table 0J Proceeding further, we can now place all 
pairs of the three basis states on a single Mn-Mn bond and form states of given total 
bond spin Sh- Then we need to calculate the overlap to all pairs of excited states that 
are connected to our bond ground-state by a single electron transfer and share the 
same bond spin Sb- Clearly, there is a reasonable number of different combinations, 
and the use of some computer algebra system is recommended. Each of the matrix 
elements 

{d n d m \H hovv \d n - 1 d m+1 ) = O a pS{S l S ] )t M , a,(3e{6,e,a 2 } (18) 

can be decomposed into an orbital part Q( a, a spin part S(SiSj) and the corresponding 
transfer amplitude t^y In Table we list these components for all d 4 d 4 ^ d 3 d 5 , 
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d 4 d 3 ^ d 3 d 4 , d 3 d A ^ d 2 d 5 and d 3 d 3 ^ d 2 d 4 excitations for a bond in z-direction, i.e. 
j = i + z. Comparing the orbital parts O a/ 3 with the rotation properties of the 9 and 
£ functions given in Table |3 we observe that the O a p can be simplified, if we change 
the orbital quantisation axis from z to x or y. For instance, the four matrix elements 
— {3, \/3, V3, l}/4 in line 6 of Tabled which connect the functions {99, 9e, e9, ee} to 
the appropriate excitation, reduce to a single matrix element —1 originating from the 
bond state e x e x with e x = — (V39 + e)/2. Hereafter we use these properties to express 
some of the interactions with the help of the rotation operators Rs, S £ {x, y, z}. 

Putting together all the above pieces we are now in the position to formulate 
the second order contributions to our effective electronic Hamiltonian. For simplicity, 
we neglect all terms which involve hopping processes between three sites, and restrict 
ourselves to those terms, where an electron hops back and forth on a single Mn-Mn 
bond, 

l^a^^ajl ffhopp |^)(^| gho PP |A/? 3 )(A/3j| 



H = — 



{ij) 



Ae(#) 



(19) 



In addition, we follow Ref. [26] by assuming C ~ 4B and introducing new Coulomb 
parameters U and Jh, which refer to a Hubbard like on-site repulsion and a Hunds 
rule coupling, respectively. The actual choice of the relation between U, Jh and A, 
B, C = AB is a matter of convention. Having in mind the undoped compounds and 
a bond with d 4 ( 5 E) basis state on each site, we may ascribe U to the minimal energy 
required for transferring one of the e g electrons and forming a low spin d 5 state on 
one of the sites. On the other hand, the excitation with a high-spin d 5 state should 
have energy U — 5Jh, where the factor 5 reminds of the energy difference originating 
from an exchange term JhSs with \S\ = 2 and |s| = 1/2. This choice leads to 

U = e[d 5 (%)d 3 (%)] - e[d 4 ( 5 E)d 4 ( 5 E)} = A + 2B + 5C , 
5J h - eidH^d 3 ^)} -e[d 5 (%)d 3 (%)] = WB + 5C, (20) 
and after some algebra for a z-bond we obtain the following contribution of second 
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H.c. 



(21) 



The first three terms reproduce the result of Ref. [26] for the undoped compounds, 
which are characterised by a competition of one ferromagnetic and a number of 
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Figure 4. Elongation patterns of the oxygen ions within the E g and A\ g vibration 
modes of a MnOe octahedron. 

anti-ferromagnetic spin exchange terms of order f?,/f7. The latter, however, are 
directly coupled to pairs of orbital projectors favouring different patterns of orbital 
ordering. Upon doping the ferromagnetic double-exchange interaction, Equation {T7JI, 
comes into play, but also anti-ferromagnetic terms of order i^/J^ gain importance. 
Obviously, the complex phase diagram and the different ordering patterns of the 
manganites are a direct consequence of effects noticeable already on the level of the 
microscopic Hamiltonian. 

5. Electron-phonon interaction 

As discussed before, most undoped or weakly doped manganites show long range 
(Jahn- Teller type) distortions from the ideal perovskite structure, and changes in 
the magnetic or orbital ordering patterns are usually accompanied by changes in the 
crystal structure. In addition, modern experimental techniques show that besides 
these cooperative static effects also dynamic and local properties of the lattice are 
important for an understanding of the manganites. To name only a few, we may 
think of x-ray-absorption fine-structure data [31] showing a clear relationship between 
magnetic order and local lattice distortion, or measurements of the atomic pair- 
correlation function [32, 33] proving the correlation of electric resistivity and local 
lattice structure near the metal-insulator transition. It is thus natural to consider 
electron lattice interactions and, in particular, polaronic effects as an important 
ingredient of a microscopic description of CMR manganites [34,35]. 

Like for the electronic, spin and orbital degrees of freedom, the local environment 
of the manganese ions is the appropriate starting point also for the modelling of the 
electron lattice interaction. In Figure 01 we show the basic vibration modes of the 
Mn06 octahedron that are even under inversion and thus susceptible for a linear 
coupling to the d orbitals. The modes Qg and Q £ , which belong to the irreducible 
representation E g , are responsible for the Jahn- Teller effect [36,37], where the system 
tries to lift the degeneracy of the two 5 E ionic basis states (see Equation IjlOjll and 
gain electronic energy by lowering the point symmetry through a distortion of the 
lattice. As above the electronic degrees of freedom are described by the fermionic 
operators d!gj £ , whereas for the phonons we have the bosonic operators bgJ E . The 
interaction between both should be linear in the bosons, bilinear in the fermions, and 
belong to the irreducible representation A\ to form a Hamiltonian that conserves cubic 
symmetry. These requirements lead to one of the classical Jahn- Teller problems, the 
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so-called E <g> e model [38, 39] 



IT 



(n E - ra e )(4 + 6 fl ) + (djd c + 4d fl )(6| + b c ) 



btbn + bib; 



(22) 



where the last term refers to the usual harmonic lattice dynamics. For a single 
ion this model possesses an additional 0(2) symmetry, which becomes manifest 
in the commutation of the Hamiltonian _ff JT with the A2 symmetric operator 
(dgd £ — d\d e )/2 — (bgb £ — b e b e ), and leads to a two- fold degeneracy of every eigenvalue. 
Whereas for a single ion this seeming contradiction to the common notion of a Jahn- 
Teller system can be lifted by higher order electron lattice interactions, for realistic 
compounds like the manganites this symmetry is already broken by the orbitally 
anisotropic hopping between neighbouring sites, Equation lfH|l . Unfortunately we 
thereby loose an advantage of this symmetry: the simple tridiagonal structure of the 
E ® e Hamiltonian matrix [38] . Compared to a symmetric hopping, the symmetry- 
breaking orbital anisotropy may also increases the susceptibility of the system to 
Jahn- Teller type ordering and polaronic effects, an issue which could be interesting 
for future studies. 

Coming back to the phonon modes of the MnOg octahedron, the symmetric 
mode Q ai should couple to a bilinear fermionic operator belonging to the same A\ 
representation. The most natural operator of this type is the electronic density, and 
we arrive at the Holstein type [40, 41] Hamiltonian, 

Ho1 = g K + n e )(bl +b ai )+" bl b ai , (23) 

which is well known from polaron physics. Since we use hole operators instead of 
electrons, we prefer to adjust the expression for the local density such that the 5 E 
states couple to Q ai but does not, n i g + n i e — 2n i g n i £ . Generalising the above 
types of local electron phonon interaction to the crystal and neglecting the weak 
dispersion of the optical phonons corresponding to Qg, Q e and Q ai , we finally arrive 
at the Hamiltonian 



H 



EP 



[( n i,e ~ n ifi){b\ft + Ke) + - 1 

i 

'^( n ifi + n i,e ~ 2n i,e n i,e)( b i, ai + KaJ 



>)(&[, 



(24) 



-V^ b- 

/ j 2,ai l,ai 



6. Discussion of the microscopic model and numerical results 

Summarising the preceding sections, the complete microscopic model is given by 

H = J2 MHl'* +5 + Hfo s ) + EP , (25) 

i,S 

which contains all essential features of the low-energy physics of the manganites. The 
first term accounts for the ferromagnetic double-exchange interaction, which describes 
itinerant e g electrons that can optimise their kinetic energy by ordering the spin 
background formed out of localised t2 g electrons. The second term couples this spin 
background to the orbital degrees of freedom of the e g electrons via Heisenberg type 
exchange interactions modulated by orbital projectors. It is thus responsible for anti- 
ferromagnetic phases, different orbital orderings as well as spin-orbital correlations. 




Figure 5. Ground-state properties of the full model, Equation 1251 . on a four 
site cluster at doping x = 0. (a) Evolution of lattice correlations with increasing 
electron phonon (EP) coupling, (b) Phase diagram of the model without EP 
coupling; the blue star and arrow mark parameters U, used in the other panels, 
(c) Spin and orbital ordering with increasing EP coupling g. (d) Spin and orbital 
ordering with increasing U and g = 0. (e) Decoupling of spins and orbitals with 
increasing EP coupling g. (f) Spin-orbital correlations at g = 0. Throughout 
J h = 0.7 eV, t = 3t w = 0.4 eV, u = u> = 0.07 eV, g = g; see [42] for more details. 




Figure 6. Ground-state properties of the full model, Equation 1251 . on a four 
site cluster at doping x = 0.25. (a) Due to polaron formation the kinetic energy 
decreases with g. (b) At strong g double-exchange is too weak to promote 
ferromagnetism, i.e. the total spin of the clusters goes to zero, (c) Bond-lengths 
(ix/y) along x and y direction show lattice distortion in the spin disordered phase, 
(d) The orbital orientation ip in the neighbourhood of a hole is orbital-polaron 
like only for mobile carriers, (e) and (f) Electron- phonon interaction causes an 
effective decoupling of spin and orbital degrees of freedom, (g) and (h) For 
increasing g lattice fluctuations show a kink near the FM-AFM transition, quite 
similar to the experimentally observed behaviour as a function of temperature [31] . 
See [42,43] for more details. 
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Figure 7. Ground-state properties of the full model, Equation 1251 . on a four site 
cluster at doping x = 0.5. (a) Similar to x = 0.25 increasing g causes a polaronic 
band narrowing, but due to charge ordering and more important AF interactions 
magnetism is decoupled from charge dynamics, (b) For small g orbital correlations 
can be characterised by complex orbital states, (c) Strong EP coupling leads to 
charge order, (d) The lattice can show different distortion patterns. See [42] for 
more details. 




Figure 8. Schematic summary of the various spin, charge, orbital and lattice 
patterns calculated within the full model, Equation (1251 ■ at weak and strong 
electron phonon interaction and increasing doping. 
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The third term, finally, causes all the polaronic effects and long range lattice distortions 
observed in some regions of the manganite phase diagram, but also affects the orbital 
ordering. 

For analytical methods the above Hamiltonian is far too complex to be understood 
in full detail, and even its numerical solution is hard. Using high performance 
computers and new optimisation methods for the lattice degrees of freedom, we 
were able to study the ground-state properties of small clusters and to address, 
in particular, short range correlations [42,43]. In Figures 03 to 03 we give a short 
overview of these results. The main features are the following: In the undoped 
compounds, e.g. LaMnC>3, only the second and third term of Equation (|2~5|) are 
active, and without electron phonon interaction the competition of the spin-orbital 
contributions depends sensitively on the values of the Coulomb and Hunds rule 
coupling U and J^, respectively. The dynamics of both subsystems is strongly 
correlated. With increasing electron phonon interaction the systems tends to develop 
static Jahn- Teller distortions, which also fix the orbital pattern and subsequently the 
spin order. Dynamic correlations of spins and orbitals are suppressed. Upon doping 
the ferromagnetic double-exchange comes into play, and only a substantial polaronic 
band narrowing due to strong electron lattice interaction can prevent ferromagnetic 
order. The orbital dynamics is coupled mainly to the charge dynamics, and is less 
affected by the spin-orbital terms (see Figure El where also the case U = oo, i.e. the 
disabling of most of the spin-orbital terms, is considered). Interestingly, changes 
in the spin order are reflected also in the lattice fluctuations, an effect observed 
as well in experiment [31]. Facilitated by the electron lattice interaction at and 
above doping x = 0.5 the susceptibility of the system to charge ordering dominates 
the low-energy behaviour and matching spin interactions lead to antiferromagnetic 
order. These features are, of course, well known from experiment [44-46]. The 
orbital correlations can show interesting patterns, involving, for instance, complex 
linear combinations \<fr) = cos(4>)\9) + sm(4>)\e) of the e g states combined to two-site 
correlations proportional to \4>}i ® \4>')j + \4>')i ® \4>)j with nonzero imaginary parts of 
4> and if>'. Similar complex correlations were also found in mean-field type treatments 
of the manganites [47-49] . 

Of course, numerical studies of small clusters are of limited value if we want to 
address long-range order or thermodynamic properties. For the former, mean-field 
studies of the complete microscopic model may give some insight and were carried 
out e.g. for the undoped compounds [26]. However, since a large number of ordering 
patterns usually differs by only tiny amounts of energy, the results of such calculations 
often depend on the underlying assumptions. The study of thermodynamic properties, 
on the other hand, requires further simplifications of the model, which will depend on 
the particular aspect of the manganite system that we are interested in. Hereafter we 
discuss some candidates of such simplified models. 

One of the first models that was studied in connexion with the manganites is the 
double-exchange Hamiltonian [29,30,50], 

H DE = ]T (aU*r 4cj + H.c.) with J2 4r°i = 2S + °i c i Vi > ( 26 ) 

which follows from the first term in Equation l|25|l by omitting the orbital degrees 
of freedom and generalising to arbitrary on-site spin S. Despite its apparently 
simple structure - bosons coupled to spinless fermions - the model is not solvable 
directly, but many of its properties can been studied within the coherent potential 
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Figure 9. Density of non-zero eigenvalues of the quantum double-exchange model 
with spins of amplitude S (dotted line: discrete spectrum; green dot-dashed line: 
running average) compared to the classical double-exchange model (bold red line). 
Left: Two fermions on a four site ring. Right: Four fermions on an eight site ring. 
See [28] for more details. 



approximation (CPA) [51-53] or, even more basic, by mean-field approximation [50]. 
The latter can be obtained by considering the matrix element for a single bond [30], 
t = t(ST + l/2)/(25 + 1), and averaging over all orientations of the total bond spin 
St within an effective mean field A = [3gfiBH* s [28,50]. The double-exchange model 
then reduces to the non-interacting fcrmion Hamiltonian 



H. 



DE,1 
off 



-t 7s [S'A]£( C ic j +H.c 
coth((5 



5 + 1 



1)A) 



25+1 



with 



coth((5 + 1/2)A) - 



(27) 



coth(A/2) 
25+1 



and for a given temperature the parameter A is chosen to minimise the free energy of 
the fermions and the spin system. 

If the spin background is less important and the focus is more on the electrons 
and their interaction with other degrees of freedom (e.g. phonons, as discussed below), 
we may completely neglect the quantum nature of the spins. The most direct way to 
find the limit 5 — > oo of Equation (|26|l is the average over spin coherent states [54] 
which, for a classical spin background parameterised by the polar angles {0k,4>k}, 
yields [55,56] 



H: 



DE 



(ij) 



t l3 c\ Cj + H.c. 



with 



(28) 



To access the quality of this approximation, in a recent work [28] we calculated the 
density of states (DOS) of the double-exchange model for small finite systems and 
compared the cases of quantum and classical spins. As Figure illustrates, at least 
in a disordered spin background the classical approximation describes the spectrum 
very well, even in the case 5 = 3/2 that is relevant for the manganites. 

In analogy to the mean-field model, Equation 127(1 , where the spin background is 
taken into account only on average, we can also average over the fermion degrees of 
freedom to obtain an effective spin Hamiltonian. For the case of classical spins this 
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Figure 10. (a) Comparison of the magnetisation M(T) of the classical double- 
exchange model, Equation 1281 . and of the effective model, Equation 1291 . 
obtained from Monte Carlo simulations [57, 58] of L 3 site clusters, (b) and (c) 
Suppression of the ferromagnetic order by electron phonon interaction; e p = g 2 /u) 
denotes the coupling strength of the underlying Holstein double-exchange model. 
See [58] for more details. 

leads to 

#eff E ' 2 = - J off \A+ ^ With Je& = tx(l-X)/V2. (29) 

(ij) 

Monte Carlo simulations [57, 58] show that the magnetisation data and critical 
temperatures of this model agree surprisingly well with those for the full 
classical double-exchange model [see Figure ITU1 fa)]. Summarising these different 
approximations we arrive at the conclusion that the double-exchange is well described 
by simple effective models, in particular, if we consider manganites with their almost 
classical spin amplitude S = 3/2. 

As discussed above, the double-exchange interaction is responsible for the 
ferromagnetic metallic phase of the manganites. However, soon after the discovery 
of the colossal magneto-resistance effect additional electron-lattice interactions were 
realised to be important for an understanding of the conductivity [34,35]. The third 
part of our complete model, £f EP , contains three different phonon modes that couple 
to the orbital degrees of freedom or to the density, but certainly not all of them are 
equally important in different regions of the manganite phase diagram. For instance, 
matching the double-exchange model in Equation l|26fl . we may again neglect orbital 
degrees of freedom and the corresponding phonon modes, which directly leads to the 
Holstein double-exchange model 

H = 2sT\ £ c ^ + H - c -) + 9 + b ^ +U E b * b > • ( 3 °) 

(ij).o i i 

Its properties have been studied with a number of methods and seem to account 
for many aspects of the metal-insulator transition and of the transport properties of 
ferromagnetic-metallic manganites with 0.25 < x < 0.5. A good overview is given e.g. 
in Ref. [59] with a focus on CPA. Analogous to the pure double-exchange Hamiltonian, 
we can also apply a number of further approximations to Equation l|30l) . e.g. classical 
spins and phonons or Lang-Firsov type variational treatments of the lattice dynamics. 
In Figure El (b) and (c) we illustrate the suppression of ferromagnetic order by the 
electron-phonon interaction. To obtain these results the spin part of Equation (|30ll is 
approximated by the effective spin model from Equation (|29|l and the electron-phonon 
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part is treated within a modified variational Lang-Firsov approximation [60,61]. The 
resulting classical model is simulated with a Cluster Monte Carlo Method [58] . Finally 
note, that also the polaron problems related to _ff EP are interesting by themselves. 
Neglecting the spin degrees of freedom, our full Hamiltonian (|25[1 reduces to a fermionic 
hopping term plus electron phonon interactions of Holstein and E®e Jahn- Teller type. 
The corresponding few-fermion problems define the Holstein polaron and E ® e Jahn- 
Teller polaron, whose properties have been studied numerically [62,63]. 

Having discussed simplifications mainly of the double-exchange (first) and polaron 
(third) part of our microscopic model, we left out the spin-orbital (second) part so 
far. However, except for the various mean-field treatments, a further simplification, 
which does not alter the delicate balance of these degrees of freedom, is less obvious. 
On the one hand, there are a number of numerical studies [64-66] which focus on the 
orbital dynamics within a fixed spin background. On the other hand, we can study 
the interplay of spins and orbitals from a purely theoretical point of view and consider 
simplified models like that derived by Kugel and Khomskii [67] 



4(S,S,)(r/ + l)(r/ + \) + (rf - i)(r/ - \) 



(31) 



where the orbital pseudo spin operators t s are related to our orbital projectors via 
R s (P e / £ ) = 1/2 ± t 5 . Recently, this model was considered in connexion to basic 
mechanisms of spin-orbital interaction, e.g. quantum disorder versus order-out-of- 
disorder [68-70]. Note also that there is a one-dimensional variant of the model, 
which is exactly solvable [71,72]. 



7. Summary 

In this contribution we focussed on a detailed introduction to the microscopic 
modelling of colossal magneto-resistance manganites, a class of materials which is 
particularly interesting due the close interplay of charge, spin, orbital and lattice 
degrees of freedom characterising its phase diagram and electronic properties. We 
supplemented the derivation by a discussion of various simplified models and their 
relation to our Hamiltonian. In addition, we reviewed some of our previous numerical 
results for the complete microscopic model and for its relatives. 
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